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Abstract 

Results from electric charge correlations studied with the Balance 
Function method in A+A collisions from 20A to 158^ GcV are pre- 
sented in two different rapidity intervals: In the mid-rapidity region 
we observe a decrease of the width of the Balance Function distri- 
bution with increasing centrality of the collision, whereas this effect 
vanishes in the forward rapidity region. 

Results from the energy dependence study in central Pb-|-Pb col- 
lisions show that the narrowing of the Balance Function expressed by 
the normalised width parameter W increases with energy towards the 
highest SPS and RHIC energies. 

Finally we compare our experimental data points with predictions 
of several models. The hadronic string models UrQMD and HIJING 
do not reproduce the observed narrowing of the Balance Function. 
However, AMPT which contains a quark-parton transport phase be- 
fore hadronization can reproduce the narrowing of the BF's width 
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with centrality. This confirms the proposed sensitivity of the Balance 
Function analysis to the time of hadronization. 



1 Introduction 

The study of the Balance Function (BF) is motivated by the prediction 
that its width may be sensitive to the creation of a deconfined phase in the 
early stage of central nucleus- nucleus collisions [HEllS]. The method is based 
on the observation that hadrons are produced locally in pairs of oppositely 
charged particle. Particles of such a pair are separated in rapidity (y) due 
to the initial momentum difference and secondary interactions with other 
particles. Particles of a pair that were created earlier are separated further 
in rapidity because of the expected large initial momentum difference and 
the long lasting rescattering phase. On the other hand, oppositely charged 
particles of a pair that were created later are correlated over a smaller interval 
A?/ in rapidity. 

As in previous publications from the NA49 experiment on the BF 
the analysis is performed for all charged particles (predominantly pions) in 
order to avoid the additional complications from the identification procedure. 
Consequently the pseudo-rapidity (rj) of charged particles is examined and 
the BF is defined as [T]: 



B{Ari) 



N+_ {At]) -N_. ( Ar/) A^_+ {At]) - N++ ( Ar/) 



+ 



A^4 



The width of the BF can be characterized by the weighted average {At]) : 



{Ar^) = j2iB.-^m)/j:B„ (2) 

1=0 1=0 

where i is the bin number of the BF histogram. 

It was suggested that the creation of a deconfined phase of quarks and 
gluons (the QGP) in the early stages of nucleus-nucleus collisions, would lead 
to delayed hadronization, i.e. a narrowing of the BF Qjj. Indeed, results on 
the BF obtained for Au+Au collisions at RHIC [3] and Pb-|-Pb interactions 
at the top SPS energy [1] observed this narrowing for central collisions in 
the mid-rapidity domain. In this paper, the study of the BF is extended to 



3 



the forward rapidity region. In addition, the energy dependence of the width 
of the BF for central Pb+Pb colhsions is investigated at lower SPS energies 
where indications for the onset of deconfinement had been obtained [6]. 

In the following two sections we will first describe in brief the NA49 
experimental setup and the data analysis. The next two sections will be 
dedicated to the new experimental results on the rapidity and the energy 
dependence of the width of the BF. We will conclude with a section on 
model comparisons and the summary. 

2 Experimental Setup 

The NA49 detector [7J is a wide acceptance hadron spectrometer for the 
study of hadron production in collisions of hadrons or heavy ions at the 
CERN SPS. The main components are four large - volume Time Projection 
Chambers (TPCs) which are capable of detecting 80% of some 1500 charged 
particles created in a central Pb+Pb collision at 158^4 GeV. 

The targets are C (561 mg/cm^), Si (1170 mg/cm^) disks and a Pb (224 
mg/cm^) foil for ion collisions and a liquid hydrogen cylinder (length 20 cm) 
for hadron interactions. They are positioned about 80 cm upstream from 
VTPC-1. 

The centrality of a collision is selected (on-line for central Pb+Pb, Si+Si 
and C+C and off-line for minimum bias Pb+Pb, Si+Si and C+C interac- 
tions) by a trigger using information from a downstream calorimeter (VCAL), 
which measures the energy Eq of the projectile spectator nucleons. 

3 Data Analysis 

3.1 Data Sets 

The data sets used in the rapidity dependence study come from p+p, 
C+C, Si+Si and Pb+Pb collisions at 158^ GeV (Table [T]) and at 40^ GeV 
(Table [2]). For Pb+Pb interactions, data with both central and minimum 
bias trigger have been analyzed in order to study the centrality dependence of 
the BF. The minimum bias data were subdivided into six different centrality 
classes [BJ according to the energy recorded by the VCAL, from class Veto 1 
(the most central collisions) to class Veto 6 (the most peripheral collisions) 
(Table [Hand [2]). For the energy dependence study, we have analyzed central 
Pb+Pb collisions throughout the whole SPS energy range (Table [3]). The 
most central Pb+Pb interactions for the highest SPS energy correspond to 
5% whereas for lower energies to 7.2% of the total geometric cross section. 
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Interaction 


Nnnil)or of events 




(A'n) 


P+P 


200K 






c+c 


45K 


15.3 


14 


Si+Si 


65K 


12.2 


37 


Pb+Pb(6) 


38K 


43.5 - ... 


45 ±4 


Pb+Pb(5) 


20K 


33.5 - 43.5 


85 ±7 


Pb+Pb(4) 


24K 


23.5 - 33.5 


128 ±8 


Pb+Pb(3) 


24K 


12.5 - 23.5 


196 ±6 


Pb+Pb(2) 


50K 


5 - 12.5 


281 ±4 


Pb+Pb(l) 


38K 


- 5 


352 ±3 



Tabic 1: Systems and centrality classes used in this analysis. Listed for 
p+p, C+C, Si+Si and six centralities of Pb+Pb collisions at 158^ GeV are 
the number of events, the cross section of selected central interactions as 
percentage of the total inelastic cross section and the mean number {Nw) of 
wounded nucleons. 

3.2 Event and Track Selection 

In order to reduce the contamination from non-target events and non- 
vertex tracks, selection criteria were imposed both at the event and the track 
level. 

Events were selected that had a proper position of the reconstructed pri- 
mary vertex. The vertex coordinate Vz along the beam axis had to fulfill 
1^ — ^^ol < where the values of the central position V^g and the range 

vary with respect to the system and the energy analyzed. In addition the 
vertex coordinates \4 and Vy perpendicular to the beam axis had to fulfill 
IK: — Kol < and — V^ol < ^Vi where the values V^^, Vy^ and Ax, Ay 
also vary accordingly. 

Selection criteria at the track level were imposed in order to reduce the 
contamination by tracks from weak decays, secondary interactions and other 
sources of non- vertex tracks. Thus, an accepted track had to have an extrap- 
olated distance of closest approach d,j. and dy of the particle at the vertex 
plane within the range \dx\ < 2.0 cm and \dy\ < 1.0 cm. In addition the 
potential number of points in the detector for the selected tracks had to be 
more than 30. To suppress double counting due to track splitting the ratio 
of the number of reconstructed points to the potential number of points was 
required to be larger than 0.5. 

Finally, the NA49 detectors provide large acceptance in momentum space; 
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Interaction 


Number of events 




{Nw) 


C+C 


76K 


65.7 


9 


Si+Si 


43K 


29.2 


32 


Pb+Pb(6) 


45K 


43.5 - 


43 ±4 


Pb+Pb(5) 


27K 


33.5 - 43.5 


93 ±7 


Pb+Pb(4) 


30K 


23.5 - 33.5 


142 ±8 


Pb+Pb(3) 


32K 


12.5 - 23.5 


210 ±6 


Pb+Pb(2) 


106K 


5 - 12.5 


290 ±4 


Pb+Pb(l) 


190K 


- 5 


351 ±3 



Table 2: Systems and centrality classes used in this analysis. Listed for C+C, 
Si+Si and six centralities of Pb+Pb collisions at 40^4 GeV are the number of 
events, the cross section of selected central interactions as percentage of the 
total inelastic cross section and the mean number {Nw) of wounded nucleons. 

however the acceptance in the azimuthal angle is not complete. To avoid 
edge effects and to facilitate comparisons with model calculations, particles 
were selected in a restricted acceptance region [9j. The boundary of this 
region is described by Eq. [SJ 

M'^) = ^ + (L)a + (3) 

where the values of the parameters A, B, C and D depend on the rapidity 
interval and on the energy and can be found in [H]. 

4 Rapidity Dependence 

In order to investigate the properties of hadronization in heavy ion colli- 
sions, as proposed by the BF methodology [T], we have studied the system size 
and centrality dependence of the width of the BF at two SPS energies (158^4 
and 40A GeV, corresponding to ^/snn = 17.3 and 8.8 GeV, respectively) 
and in two different pseudo-rapidity intervals. The selected pseudo-rapidity 
regions for ^/s = 17.3 GeV are 2.5 < r] < 3.9 named as mid-rapidity region 
and 4.0 < f] < 5.4 named as forward rapidity region. The corresponding 
pseudo-rapidity regions for y/s = 8.8 GeV are 1.8 < r] < 3.2 (mid-rapidity 
region) and 3.3 < rj < 4.7 (forward rapidity region). Fig. [T] shows the re- 
sulting BF distributions for the two different energies (upper and lower plots 
respectively) as a function of the pseudo-rapidity difference Ar/ for central 
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Figure 1: The BF as a function of the pseudo-rapidity difference Arj for cen- 
tral (centrahty class 1 of Tables [Hand [2]- left plots) and peripheral (centrality 
class 6 of Tables [Hand [2]- right plots) Pb+Pb colhsions at ^/sNN = 17.3 GeV 
(upper plots) and ^/snn = 8.8 GeV (lower plots) for the two rapidity regions. 
The curves drawn are Gaussian fits used to guide the eye. 
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Figure 2: The system size and centrality dependence of the measured width 
{At]) of the Balance Function for charged particles at ^/sNN = 17.3 GeV (top 
plots) and ^/snn = 8.8 GeV (bottom plots) plotted as a function of the mean 
number of wounded nucleons (Ny^) for two different rapidity regions (Table 
HI): the mid-rapidity region (left plots) and the forward rapidity region (right 
plots). 
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y/SNN [^eV] 


Ef^arr. [A GcV] 


O'centJ 0"me/.[%] 


{Nw) 


6.3 


20 


7 


349 ± 1 ± 5 


7.6 


30 


7 


349 ± 1 ± 5 


8.8 


40 


5 


361 ±1 ±5 


12.3 


80 


7 


349 ± 1 ± 5 


17.3 


158 


5 


362 ±8 



Table 3: The different energies for Pb+Pb interactions used to perform tlie 
energy scan. Tlie cms energy of tlie system as well as the different beam 
energies, the cross section of selected central interactions as percentage of 
the total inelastic cross section as well as the number of wounded nucleons 
are listed. The third number in the last column refers to the statistical error 
whereas the second to the systematic. 

(left plots) and peripheral (right plots) Pb+Pb collisions. 

The top plots of fig. [2] shows the width of the BF distributions for ex- 
perimental data as a function of the mean number of wounded nucleons, for 
the two pseudo-rapidity regions (mid-rapidity region-left plot, forward rapid- 
ity region-right plot) for the y/s^N = 17.3 GeV case. The figure indicates 
that there is a clear centrality dependence for real data in the mid-rapidity 
region while this dependence vanishes in the forward region. In order to fur- 
ther investigate this effect, the Ultra-relativistic Quantum Molecular Dynam- 
ics model (UrQMD) [lOj, which is a microscopic model of (ultra) relativistic 
heavy ion collisions in the energy range from Bevalac and SIS up to AGS, 
SPS and RHIC, was used. The results from the study of this conventional 
hadronic model as well as the widths obtained from shufflecfl data [HI E] are 
also plotted in fig. [21 No system size and centrality dependence, regardless 
of the rapidity region, is seen for either UrQMD or shuffled events. The same 
behavior is observed for ^/s^n = 8.8 GeV in the bottom plots of fig. [21 

One should also note that there is no significant difference in the actual 
values of the widths between the two rapidity regions at both energies for 
the small systems, such as p+p up to Si+Si. In other words, the difference 
in the absolute values appears only for the most central Pb+Pb collisions at 
both energies. 

"'^In shuffled events the particle (pseudo)rapidities in the event are randomly reassigned 
while the particle charges are unchanged. This procedure destroys (pseudo)rapidity cor- 
relations but retains global charge conservation effects. 
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5 Energy Dependence 



y/SNN[(^eV\ 


rj interval (mid-rapidity) 


7] (forward rapidity) 


6.3 


1.6 < r/ < 3.0 


3.1 < 7] < 4.5 


7.6 


1.7 < r/ < 3.1 


3.2<r]< 4.6 


8.8 


1.8<r]< 3.2 


3.3 < ?7 < 4.7 


12.3 


2.2 < 7] < 3.6 


3.7 < < 5.1 


17.3 


2.5 < 7] < 3.9 


4.0 < 7] < 5.4 



Table 4: The table lists the definition of the mid-rapidity as well as the 
forward rapidity interval for the different SPS energies. 

The energy dependence of the BF was studied for central Pb+Pb collisions 
for which NA49 data is available throughout the SPS energy range (lOOK 
events were analyzed for each energy). These data samples passed once 
again through the shuffling mechanism to provide an estimate of the highest 
value of the width for each energy under the constraint of global charge 
conservation. The pseudo-rapidity interval analyzed was limited to a range 
of 1.4 units and was located around mid-rapidity (Table Hj). In order to 
compare the decrease of the width for the different energies, we introduced 
the normalized parameter W which is defined by the following equation: 



^ 100 ■ {{/\vi) shuffled - {^V)data) 
shuffled 

By using the measure W, we try to quantify the change of the BF width 
due to effects other than those of global charge conservation. Moreover, the 
slight changes of acceptance with energy should approximately cancel. Fig. 
O shows as filled squares the BF distributions as a function of the pseudo- 
rapidity difference Arj for central Pb-fPb collisions at ^/sNN = 12.3 GeV 
(left plot), y/SNN = 7.6 GeV (middle plot) and ^/snn = 6.3 GeV (right 
plot). Fig. [1] shows the results for the other two energies ^/snn = 17.3 GeV 
and y/SNN = 8.8 GeV. 

The left plots of fig. H] shows the dependence of the normalized W param- 
eter on y/SNN- For the data (red squares in Fig. Hj), we note an indication 
for an increase with energy. In contrast, the results from the analysis of 
UrQMD generated events (blue boxes in Fig. Hj), show no significant energy 
dependence of the measure W. 
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Figure 3: The BF as a function of the pseudo-rapidity difference Arj for 
central Pb+Pb colhsions at y/s^N = 12.3 GeV (left plot), ^/snn = 7.6 GeV 
(middle plot) and ^/sNN = 6.3 GeV (right plot) for the two rapidity regions. 
The curves drawn are Gaussian fits used to guide the eye. 
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Figure 4: The dependence of the normalized parameter W on the ^/JNN for 
central Pb+Pb collisions in the SPS energy range after applying the NA49 
acceptance filter (left plot) and for central Au+Au collisions at RHIC (right 
plot). 
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The right plot of fig. H] shows the dependence of the normalized W pa- 
rameter on ^/sNN for the energy range of SPS together with the RHIC result 
at y/SNN = 130 GeV [5j|. The corresponding RHIC point tends to be even 
higher than the last SPS one, indicating a continued rise of the W parameter 
towards RHIC leaving an open question for the LHC energies. 

Finally, we study the energy dependence of the scaled parameter W of 
the BF's width in the forward rapidity region (Table H]). The corresponding 
BF distributions for central Pb+Pb collisions are shown as open squares in 
Fig. [Sat = 12.3 GeV (left plot), ^s]^ = 7.6 GeV (middle plot), 

^/sNN = 6.3 GeV (right plot) and in Fig. [T] at ^/sNN = 8.8 GeV (bottom 
left) and ^/sNN = l7.3GeV (top left). Fig. E] summarizes the main result 
from this study: The normalized parameter W shows no sign of an energy 
dependence in the forward rapidity regions. 

6 Model comparison 

In an attempt to interpret the results obtained for the centrality depen- 
dence of the BF width, several microscopic models were analyzed such as 
HIJING dH, UrQMD [lO] and AMPT The first is based on the exci- 
tation of strings and their subsequent hadronization according to the LUND 
model [13]. The second, is also based on string excitation with the addi- 
tion of interactions between the produced particles after the hadronization 
(rescattering) . The last is a multi-phase transport model which contains 
a quark-parton transport phase before hadronization. In AMPT, the initial 
spatial and momentum distributions of partons and string excitations are ob- 
tained from HIJING. The parton cascade follows the ZPC model When 
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the interactions of the partons cease, they are recombined with their parent 
strings to form hadrons according to the LUND fragmentation mechanism 
[T3] . Then the scatterings of the produced hadrons are described by the ART 
model [T5] . 

Figure El in which the measured width {At]) of the BF for charged par- 
ticles at y/S]\fN = 17.3 GeV is plotted as a function of the mean number 
of wounded nucleons (N^), shows the main results of our comparison. The 
clear centrality dependence of the order of 17 ± 3% seen for experimental 
data, cannot be described by the microscopic models HIJING and UrQMD. 
Results from both models suggest no centrality dependence of the BF width. 
On the other hand, the AMPT model, which incorporates the time evolution 
of partons before their hadronization, is predicting a centrality dependence 
of the BF similar to what is seen in the data. 



Figure 6: The centrality dependence of 
the measured width {Ar]) of the Bal- 
ance Function for charged particles at 
^/sNN = 17.3 GeV as a function of 
the mean number of wounded nucle- 
ons (Ny^). The plot shows the ex- 
perimental data along with the data 
points from different microscopic mod- 
els. The interval analyzed can be seen 
in Table H 



More comprehensive model comparisons are in progress and will be ad- 
dressed in a future publication. 

7 Summary 

In this paper, we have presented new results on the BF obtained by the 
NA49 collaboration. We have analyzed data from different colliding systems 
and centrality classes for two different energies in two different rapidity re- 
gions. At mid-rapidity the experimental data reveal an interesting decrease 
of the BF width with increasing centrality of the collision for both energies. 
This effect can not be reproduced either by the UrQMD model or by the 
shuffling mechanism. In the forward rapidity region, on the other hand, no 
centrality dependence of the BF width is observed. 
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The results from the energy scan show an indication for an energy de- 
pendence which is not apparent in the UrQMD model. The narrowing of the 
Balance Function, expressed by the normalised width parameter W, increases 
from the lowest SPS energy to the highest SPS and RHIC energies. Again 
this rise is limited to the mid-rapidity regions. 

Finally, the comparison with different microscopic models, showed that 
the narrowing of the BF's width with centrality can only be reproduced by 
models that contain a quark-parton evolution phase before hadronization. 
The latter might be a first indication that indeed the BF is related to the 
time of hadronization. 

This paper is dedicated to the memory of Angelos Petridis. 
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